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Abstract 

We study the finite size spectrum of integrable quantum chains of interacting non-Abelian anyons 
constructed using the Drinfeld double of the dihedral group D3 . The gapless low energy modes are 
identified as the direct product of two conformal field theories which can be decomposed according 
to the residual symmetries of the chains subject to periodic boundary conditions. 
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Ground state degeneracies in a many particle system may be an indication for the pres- 
ence of topological order without a local order parameter. Well known examples for such 
topological quantum liquids are the fractional Hall states [lj]. Another class of systems 
where transitions between quantum phases driven by topology have been proposed are two- 



dimensional frustrated quantum magnets [2 



Possible realizations for these spin-liquid 



states are certain Iridium compounds [5]. As a first step towards the characterization of 
the different phases in a given model the quasi-particles of the corresponding theory need 
to be identified. In a (2 + 1) topological quantum liquid these collective excitations can 
be described as defects in a planar gauge theory in a broken symmetry phase with finite 
residual gauge group H |6|: the quasi-particles in these systems are irreducible representa- 
tions (irreps) of the Drinfeld double D(H) labeled by their flux, i.e. an element of h G H, 
and their topological charge determined by the transformation properties under the resid- 
ual global symmetry commuting with the flux h. States with several quasi-particles can 
be manipulated by two operations: (1) fusion is determined by the decomposition of prod- 
uct states into irreps of the Drinfeld double. (2) Interchange of two constituents in this 
state may reveal non-trivial anyonic statistics corresponding to a representation of the braid 
group. Hence, representations of Drinfeld doubles allow for the description of particles with 
anyonic statistics: apart from the appearance of a phase factor, as in the case of Abelian 
anyons, braiding may correspond to a unitary rotation of the wave function in a degenerate 
manifold. The fact that these non- Abelian anyons are protected by their topological charge 
makes them potentially interesting as resources for quantum computation [7|, |8|. 

This picture provides a framework for the investigation of topological phases and quan- 
tum phase transitions : Microscopic lattice models satisfying these constraints have been 
obtained for the non-Abelian deg rees of freedom in su(2)k Chern Simons theories, e.g. Ising 



4| and Fibonacci anyons [10|, LL3J. Local interactions between these objects can be varied to 
favour certain fusion channels. This allows to explore the phase diagram of these systems 
and to study their critical properties near quantum phase transitions. This approach works 
particularly well for quasi one- dimensional anyonic models such as chains or ladders where 
powerful numerical and analytical methods are available: at the quantum phase transition 
the low energy effective theory of these systems is expected to be a conformal field the- 
ory (CFT) and the universality class is determined by the central charge of the underlying 
Virasoro algebra. At the same time anyonic chains can be seen as realizations for the in- 



terface between phases with different topological order 12H15I] and the CFT determines the 
properties of gapless edge modes propagating along these interfaces. 

In this letter we study the critical behaviour of a system of interacting anyons within an 
integrable quantum chain model which is constructed directly using the algebraic structure 
of the gauge theory with the dihedral group of order 6 as its gauge group, i.e. the Drinfeld 
double D(D 3 ). The general representations of the Drinfeld doubles of finite group algebras 



are well known 



16]: for D(D 3 ) one has two one-dimensional irreducible representations rtf, 



four two-dimensional ones TT^' b \ (a, 6) = (0,1) or (1,6) with 6 



0,1,2, and two three- 



dimensional ones irf. This leads to an increased number of fusion channels in D(D 3 ) models 



17 in the context of 2D lattice 



as compared to the su(2)k anyon chains, see e.g. Refs. Q and 
models. On the other hand D{D$) forms a quasi-triangular Hopf algebra which provides it 
with a natural tensor product structure: the so-called quantum dimensions of these anyons 
are integers. Therefore we can consider a chain of length L with a D(D 3 ) anyon with three 
internal states represented by a copy of n 3 on each site. This spin chain is a particular 



limit of the spin 1 Fateev-Zamolodchikov model [18] and differs from the corresponding 
anyonic chain in the boundary conditions. For periodic boundary conditions the resulting 
integrable model is given by a one parametric hamiltonian with interactions between spins 
on neighbouring sites [19] 

L 

H e = J2 cos 9 h$ +1 + sin 9 hf} +1 . (1) 



In terms of the projection operators on irreps [20|| appearing in the tensor product 7r^~ <g> 7r 3 
the local hamiltonian can be expressed as 
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The operator = Ph^P = (h^)* can be obtained either by a permutation P of the spins 
on the neighbouring sites or by complex conjugation: with the exception of p^' 1 ** = (p^ 2 ^ 
the projection operators are invariant under these operations. Depending on the parameter 
9 the interactions in (JT]) favour different channels for the fusion of spins on neighbouring 
sites, e.g. the vacuum channel irf for tt < 9 < 37r/2. 

By construction the local hamiltonians have the full D(D$) symmetry. This symmetry, 
however, is broken by imposing periodic boundary conditions in the global hamiltonian ([[]), 



see Ref. 
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based on the representation theory of D(D 3 ) the tensor product (vr ; ^) ,x,i can 



be decomposed into a sum of one- and two-dimensional irreps for even L (as in Eq. (j2J) for 
L — 2). For the periodic chain ([T]) only a partial decomposition of the Hilbert space is 
possible. Below we shall use this fact to assign states to (sums of) irreps n± © 7rj~, tt^' 1 ^, 
tt^'°\ or ^ ©7T2 1 ' 2 ''. For odd L only the three-dimensional irreps 7r^ appear in the tensor 
product (vr^)® L . In the spectrum of the periodic chain these irreps cannot be separated but 
we find a connection between the eigenstates and their transformation properties under the 
residual _D 3 symmetry, namely the action of the rotation a. 

For the analysis of the spectrum of the quantum chain (pQ) we use the property that its 
hamiltonian can be split into two commuting ones with identical spectra: as a consequence 
of the decomposition (J2J) the local hamiltonians h^ 1 ' 2 ' commute as do the global ones (JT]), 
[He, He*} = 0. Furthermore, H = £\ hf^ +1 and H n / 2 = Yli hfi+i are related by a spatial 
inversion and have the same eigenvalues albeit with opposite momentum. These properties 
can be traced back to the underlying two parameter transfer matrix of this model and its 



symmetries 
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2lj . The eigenvalues of Ho = —H w can be parametrized by L complex 
rapidities Xj (u> = exp(2m/3)) 

1 — iuj e, x i 



E(X = {xA) =iV \ — coL (3) 

j 

solving the Bethe equations (j = 1, . . . , L): 

{ ' V l-iwe^ / J-l e x k _ ue x j ■ \ > 

v 7 k=l 

The spectrum of the quantum chain (pQ) is given by pairs of solutions to these equations 
corresponding to energies Ei a> p){0) = cos 6 E(X a ) + sin# E(Xp) and, similarly, momenta 
p = p(X a ) —p(Xp), provided that the combination (a, (3) satisfies the pairing rules discussed 
below. 

In the thermodynamic limit j_L ^oo, all finite solutions of Eqs. 01]) can be grouped into 
three types of so-called strings 



19|, |22J: it-strings correspond to solutions Im(xj) = 0, rr and 
2-strings are complex conjugate pairs of rapidities Xj t ± = Xj ± i2ir/3 with real center xj. 
For finite chains this classification continues to work very well for the ground states and low 
energy excitations of ±Hq. At higher energies the 2-strings become deformed, i.e. have an 
imaginary part different from ±27r/3. We have numerically diagonalized the transfer matrix 
for chains of up to L = 10 sites and found that (taking into account this deformation) all 



eigenstates can be classified this way: denoting the number of string solutions by N±, N 2 
we find iV + + N_ + 2N 2 = L — n +OQ — n_ 00 where n ±OC) e {0, 1} is the number of Bethe 
roots at x — ±00. We also find that there exist only 4 ■ 3 L / 2_1 (2 ■ 3^ L ~^' 2 ) different root 
configurations solving fll]) for even (odd) length chains. This implies that the spectrum Hg 
displays massive degeneracies (exponential in L) arising from level crossings when 6 is a 
multiple of ~. For generic values of 6 the degeneracies of the level E( a> p)(9) are lifted up to 
a remaining 'pairing multiplicity' of its components a and /3. Due to the pairing mechanism 
these components can be discussed separately. For the analysis of the low energy spectrum 
of (JTJ this amounts to the identification of the ground state and low lying excitations of 
Ha = —Ho and H . 

The ground state energy of H n has been computed in Ref. [3 for even length lattices its 
root configuration is given by a distribution of N 2 = L/2 2-strings. In the thermodynamic 
limit density functions for these strings can be introduced allowing the energy density 
to be computed. Here we have extended the analysis of the Bethe equations to excitations 
close to this state: in the corresponding configurations one or more of the 2-strings are 
replaced by ±-strings and/or Bethe roots at ±00. The spectrum of these excitations has 
a linear dispersion with Fermi velocity Vp allowing identification of the CFT for the low 



energy modes from the finite size scaling behaviour of the energies for large but finite L [231 ]: 
at the critical point the ground state energy of a 1 + 1 dimensional quantum system scales 
as Eq{L) — Leoo = —(nvp/QL) c where c is the universal central charge of the underlying 
Virasoro algebra. From the energy and momentum of low lying excitations in the finite 

system 

E(L) -E (L) = ^^(X + n + n) , p{L) - Poo = ^ (s + n - n) (5) 

the scaling dimensions X = h + h and conformal spins s = h — h of the primary fields in 
the theory can be determined (n, n are non-negative integers). 

For H^ the product Vp^c had been determined to be 12/5 previously Using Bethe 
ansatz methods and comparing the observed structure of the low energy spectrum with (jSJ) 
we find the Fermi velocity for this sector to be vp^ = 3. Hence the central charge of the 
effective field theory for the low energy degrees of freedom in H n is c = 4/5: this sector of the 
model is in the universality class of the minimal model Airs&\, the conformal weights h, h of 
the primary fields can take the rational values from the Kac table h pq = ((6p — 5g) 2 — 1)/120, 
1 < q < p < 5. The operator content of a given realization of the CFT is constrained further 



TABLE I. Scaling dimensions X n extrapolated from the finite size behaviour of the ground state and 
low energy excitations of H n for even L. (h, h) are the predictions from the .M^ g) minimal model. 
We have also indicated the D{D^) sector in which the state appears and its pairing multiplicity. 
The operator content of the sector tt^'^ is obtained from that of tt^' 1 ^ by interchanging h and h. 



D(Ds) X^ um (h, h) spin pairing mult, 

vr+evrf 0.000000(1) (0,0) 1 

0.801(3) (§,§) 1 

ttJ^ 0.4668(2) (if) -I 2 

0.666666(1) (|,0) § 2 

vrj, 1 ' 1 ) © 4 1,2) 0.13334(6) (is-is) 4 

1.33333(3) (§,§) 4 



TABLE II. As Table U for odd L. Symmetry is classified by the action of the D% 


rotation a. 


a 


jnum 


(h,h) 


spin 


pairing mult. 


1 


0.125000(5) 
0.42502(2) 
0.92490(6) 
1.625000(1) 


(o,|) 

\ 5 ' 40 > 

(2 21N 
^ 5 ' 40 1 

(o,f) 


i 

8 

3 
8 

1 

8 

13 

8 


1 
1 
1 
1 


0J, UJ~ l 


0.091665(2) 

0.59168(7) 

0.791667(1) 


fi -M 

V 15 ' 40 > 

(X 21) 

U5 ' 40 > 

(2 1) 
^3' 8' 


1 

2d 

11 

24 

13 

24 


2 
2 
2 



by modular invariance of the partition function, locality of the physical fields and boundary 
conditions 



24 



25|. 



In Tables [H [TT] we present numerical results for the scaling dimensions identified in the 
excitation spectrum of H n together with the conformal predictions. For the states given we 
have solved the Bethe equations up to a minimum of 40 sites, although in general over 100 
sites were considered when possible. Also listed are the residual symmetry sectors in which 
the levels appear. For L even the low energy spectrum of H n for a given symmetry coincides 



with that of the 3-state Potts model subject to cyclic boundary conditions with fixed Z z 



charge 



24 



26| containing spin | parafermions in the sectors rr^' 1 ^ and tc^' ^. For L odd 
'twist operators' with conformal weights h pq with q even appear in the (anti)holomorphic 
sector. The ground state is not invariant under a D% rotation. 

Note that in spite of the appearance of fields Q(h,h) with conformal spin s = h — 
h ^ Z/2 physical operators in the theory are local: they are direct products of primary 
fields h a ) (&(hp, hp)) in one-to-one correspondence to the energy levels Ei a ,p) of ([I]). 

Momentum and spin of a physical state are given by the difference between that of its two 
components. From our numerical analysis of the spectrum we find that only operators in 
the same sector with respect to the residual Z)(D3)-symmetry (the action of the rotation a) 
pair for L even (odd). With this rule the total spin of a physical field is either integer or 
half-integer. For — 7r < 9 < —tt/2 the low energy states are obtained by pairing of the states 
listed in Tables H| HT1 in the ir^'' 1 ' > sector of the model for even L the physical fields allowed by 
the pairing rules carry spin or 1, for odd L fields with spin 0, |, 1 are possible. Note that 
the observed number of different root configurations together with the pairing multiplicity 
indicated in the tables yields the total number of 3 L states of the quantum chain. 

The ground state of +H is given by a solution of the Bethe equations @ with N + = L/A 



+- strings and N_ = 3L/4 — strings [19j. It is realized for lattices of length L = (mod 4) 
and its energy scales as E (L) — Le^ = —3tt/12L [19]. The root configurations for the lowest 
excitations differ from this one by the replacement of one or two of the strings by roots at 
±oo. For excitations at higher energies 2-strings have to be taken into account. Following 



Refs. 



27H30I] the lowest finite size energy gaps are found to be 



LAE(AN±, Q±) J_ )2 _ AN+AN _ + {AN _ )2) 

271 4 o (6) 

+ l {(Q + ) 2 + Q+Q- + (Q-) 2 ) ■ 

For the lattice model the numbers AN± are related to the change in the number of ±- 
strings as compared to the ground state, i.e. take values =fL/4 (mod 1). They are further 
constrained by the fact that the total number of Bethe roots has to be L. Q± can take values 
Q± = — AN± + | (n +00 — n_oo) (mod 1). We can determine the Fermi velocity of low lying 
excitations in this sector as before finding vf,o = 3/2. Therefore the effective field theory 
for this part of the spectrum is a CFT with central charge c = 1. The field content of the 
theory is obtained from the finite size spectrum ([H]) subject to the constraints mentioned. It 



TABLE III. Scaling dimensions Xq extrapolated from the finite size behaviour of the ground state 
and low energy excitations of Hq for L = (mod 4) (the error of the extrapolation is smaller than 
the last displayed digit). (h,h) are the predictions from the Z4 parafermionic CFT. For the other 
columns, see Table HI 



D(D 3 ) 


jnum 


(h,h) 


spin 


pairing mult. 


Tit ffiTTf 


0.000000 


(0,0) 





1 




0.333332 


(o.s) 


1 
3 


2 




0.166667 


(i -M 

V12' 12 J 





4 




0.666667 







4 




TABLE IV. As Table HI] but for L = 


2 (mod 4). 




D(D 3 ) 


ynum 


(h~h) 


spin 


pairing mult. 


Tit ffiTTf 


0.750000 


(0,|) x 2,(|,0) x 2 


=c 4 


1 


^2 


0.083333 


(o,£) 


1 

12 


2 




1.083333 


U' 3' 


5 

12 


2 



TrpW^ 0.416667 ±1 4 



TABLE V. As Table IIIII for L odd. Symmetry is classified by the action of the D3 rotation a. 



a 



ynum 



(h,h) 



spin 



pairing mult. 



0.062500 
0.562500 
0.812500 



(ll,0) 

(*.o) 

U6' 4 J 



16 

_9_ 

16 

_ 11 

16 



0.145833 
0.395833 
0.645833 



(i -M 

vi6' 12/ 



16' 



V16' 12 J 



_ _1_ 

48 

_ 13 
48 

23 
18 



can be identified with that of a Z4 parafermionic theory [31|, [32J, see Tables [LTl] IIVI and IVl 
In particular, the finite size gap of the lowest states for i = L (mod 4) 7^ is determined 
by an (anti-)chiral ^=4 spin field with conformal weight ht = i(k — £)/(2k(k + 2)). 



As in the spectrum of we find states with fractional conformal spin. The physical 



8 



fields obtained after application of the pairing rules discussed above, however, have integer 
spin for < 9 < n/2 and L = (mod 4). For lattices of length L = 2 (mod 4) the spins 
can take integer or half-integer values while we find parafermions with quarter spin in 
the spectrum of chains with L odd. For 9 G (tt/2,tt) or (— ir/2,0) the operators from the 
minimal model .Mrs e) and the Z4 CFT are paired. In this regime parafermionic fields with 
quarter spin are present for L = 2 (mod 4). 

Both conformal field theories appearing in the low energy sector of the model are con- 
nected with the si (2) affine algebra by coset constructions. This allows for the interpretation 
of the quantum chain (pQ) as a description for the edge of a topological fluid nucleating within 
a surrounding non-Abelian liquid, similar to the interfaces between different quantum Hall 



states discussed recently [l2Ml5|] . As in these models the quantum critical point is pro- 
tected by topological symmetries against local perturbations. This becomes manifest in the 
small corrections to scaling fl5]) due to deviations of the quantum chain ([1]) from the CFT 
fixed point hamiltonian: the subleading L-dependence of the ground state energies indicate 
that these deviations are generated by the presence of an irrelevant operator with scaling 
dimension X = 3.82(3) in H n (X = 4.00(3) in H ). 

In summary, the exact solution of the model ([T]) has allowed to identify of the low energy 
effective theory for the critical phases in the present model in terms of the direct product of 
two CFTs. Thereby the model provides examples for interfaces between different topological 
quantum liquids supporting pairs of gapless modes. The observed pairing of the CFTs 
indicates the possibility of more general combinations of Virasoro characters in the low 



energy spectrum than those appearing in the off-diagona 
model or the Z4 parafermionic model alone 24 



modular invariants for the minimal 



25. 



321 ] . To facilitate the explicit 



construction of the new invariants appearing in the partition function of the quantum chain 



([T]) and the related D(D n ) models |2l|] the spectral analysis can be extended to different 
boundary conditions. In the context of two dimensional systems showing quantum phase 
transitions between phases with different topological order, studies of these more general 
models will be useful to obtain a CFT description of boundaries supporting several edge 
modes and their relation to the non-Abelian bulk degrees of freedom. 

Finally, we note that for the integrable cases of braided boundaries or free ends [l^ the 
full D(D n ) symmetry is restored and the models are equivalent to those obtained in the 
anyon formulation based on a fusion path along the chain. This equivalence can be used 



to embed concepts such as topological symmetry into the framework of integrable systems 
allowing for a different perspective to study possible instabilities of anyonic systems against 
local perturbations. 
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